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ABSTRACT: Characterization of the frequency response of coherent radiometric receivers is a key 
element in estimating the flux of astrophysical emissions, since the measured signal depends on the 
convolution of the source spectral emission with the instrument band shape. 
Laboratory Radio Frequency (RF) measurements of the instrument bandpass often require complex 
test setups and are subject to a number of systematic effects driven by thermal issues and impedance 
matching, particularly if cryogenic operation is involved. 

In this paper we present an approach to modeling radiometers bandpasses by integrating simu- 
lations and RF measurements of individual components. This method is based on QUCS (Quasi 
Universal Circuit Simulator), an open-source circuit simulator, which gives the flexibility of choos- 
ing among the available devices, implementing new analytical software models or using measured 
S-parameters. Therefore an independent estimate of the instrument bandpass is achieved using 
standard individual component measurements and validated analytical simulations. 
In order to automate the process of preparing input data, running simulations and exporting results 
we developed the Python package python-qucs and released it under GNU Public License . 
We discuss, as working cases, bandpass response modeling of the COFE and PLANCK Low Fre- 
quency Instrument (LFI) radiometers and compare results obtained with QUCS and with a com- 
mercial circuit simulator software. The main purpose of bandpass modeling in COFE is to optimize 
component matching, while in LFI they represent the best estimation of frequency response, since 
end-to-end measurements were strongly affected by systematic effects. 
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1. Introduction 

An accurate characterization of the bandpass response of microwave radiometers is fundamental in 
the interpretation of a wide range of astrophysical observations at millimeter wavelengths. In par- 
ticular, bandpass uncertainties may significantly impact Cosmic Microwave Background (CMB) 
experiments, since the combined emission of CMB and the Milky Way galaxy can induce system- 
atic effects both in temperature and polarization measurements. 

End-to-end measurements of the bandpasses, however, often require complex test setups that 
might be affected by systematic effects like thermal issues in the cryogenic setup or impedance 
mismatch in injecting the test signal. 

Often, a better insight of the radiometer frequency response can be achieved by building a 
bandpass model based on the S-parameters []]]] of individual components, obtained either by mea- 
surement or by analytical simulations. Such a model can give a robust independent estimate of the 
integrated radiometer bandpass, as it is based on hardware data acquired with much simpler and 
standard setups, where systematics are under control. 

The comparison of the bandpass model to end-to-end measurement or to sub-assemblies is 
valuable to identify possible unexpected interactions between the components after integration. 
Such interactions might then be modelled and included in the bandpass model. 
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1.1 Radiometers 



State-of-the-art microwave detectors in the 1-90 GHz frequency range are mainly based on coherent 
radiometers, that can achieve a sensitivity below 250 pLK^fs (at 30 GHz) thanks to the performance 
of cryogenically cooled Low Noise Amplifiers based on High Electron Mobility Transistors |@] . 

One of the most important astrophysical applications of microwave radiometers is in the ob- 
servations of the CMB. Since its first discovery in 1964[Q], an uninterrupted chain of radiometric 
observations have been carried out to measure the CMB spectral shape, anisotropy and polariza- 
tion properties with increasing precision. The most notable applications were in space missions, 
starting from COBE-DMR (COsmic Background Explorer-Differential Microwave Radiometer) 
[@] launched in 1989, with 3 channels at 31.5, 53 and 90 GHz, WMAP (Wilkinson Microwave 
Anisotropy Probe, launched in 2000) [§] with 5 channels at 23, 33, 41, 61 and 94 GHz, and the 
PLANCK Low Frequency Instrument [@], launched in 2009, with 3 radiometric channels at 30, 44, 
70 GHz. 



Feed horn 



Amplification stage Bandpass filter Diode detector 







Figure 1. Schematic of a total power radiometer 



Radiometer sensitivity scales as: 



Av 



where Av is the effective bandwidth, defined as: 



Av 



(fg(v)dv) 2 
Jg(v) 2 dv 



(1.1) 



(1.2) 



where g(v) is the frequency response. Therefore, in order to optimize the sensitivity, it is favorable 
a wide and flat band; with current technology it is reasonable to achieve typically an effective 
bandwidth of about 20% of the detection frequency, e.g. 6 GHz band at 30 GHz. 

Bandpasses do not directly impact CMB temperature anisotropy measurements, because cali- 
bration is usually performed on the CMB dipole signal which has the same black body spectrum of 
the anisotropics. However, galactic emissions have different spectral shape that will be weighted 
differently by the in-band instrument response, so that a poorly known bandpass can introduce a 
systematic error in distinguishing the CMB from the galactic emission. Furthermore, CMB po- 
larization experiments usually require differencing between independent radiometer channels with 
polarization orientations at 90° relative to each other as projected on the sky. Thus, a bandpass 
mismatch between orthogonally oriented radiometers produces a spurious differential signal that 
mixes with the sky emission, see [J7[|. 

The RF power output of the radiometer RF amplification stages is integrated into a DC volt- 
age by the detector diode, so bandpasses can only be measured with dedicated tests based on a 
monochromatic signal sweeping through the band. End-to-end frequency response tests require 
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complex setups due to cryogenics and signal injection into the horn. On the other hand, the fre- 
quency response of individual radiometer components (orthomode transducers, low noise ampli- 
fiers, waveguide components, etc.) can be tested with good accuracy and repeatability, thanks to 
the much simpler and standard setups required. It is therefore of high interest to develop methods 
to combine the measured (or simulated) information from the individual components in a software 
model of the instrument yielding a synthesised overall bandpass. 

Moreover, the new generation of CMB experiments will require focal plane arrays with hun- 
dreds or thousands of receivers in order to achieve the required sub-jitA' sensitivity per-pixel. In 
this scenario, bandpass characterization based on combination of single component data may turn 
out to be the only viable method to meet the scientific requirements. 

In this paper we present the implementation of a bandpass modeling system based on QUCS 
(Quasi Universal Circuit Simulator) 1 , a Free Software circuit simulator. The use of QUCS ensures 
the flexibility of being able to mix, whenever needed, data from laboratory measurements with 
analytically simulated components. 

Running simulations for several receivers with an interactive interface is error prone and heav- 
ily time consuming, for this reason we have developed and released under GPL 2 python-qucs 3 , 
a Python package that allows automated QUCS simulation and data exporting. 

A previous version of the bandpass model for Planck-LFI radiometers, where this work 
originated, was implemented in the ADS (Advanced Design System) software package by Agilent 
as described in [Q, gj. The use of ADS suffered several limitations: first of all, as any closed 
source software, it does not allow the customization and implementation of new features, which is 
very important for non-standard scientific applications, then the high license costs, troublesome for 
small experiments, and finally not having an easy method for scripting batch simulation runs with 
general puipose scripting languages like Python or Perl. 

In section we present a general overview on the bandpass model, in section ^3] we present 
the implementation of the rectangular waveguide component, then we show two scenarios of ap- 
plication on the COFE radiometers (section 0) and on the PLANCK radiometers (section 

2. QUCS bandpass model 

In this section we give an overview of the Quite Universal Circuit Simulator (QUCS) bandpass 
model, starting from an introduction to QUCS itself, focusing on how it is used to characterize 
bandpasses, followed by details about the rectangular waveguide model we implemented and the 
supporting Python tools used to automate the simulation and data export process. 

2.1 QUCS 

QUCS, jrO|], is an open source electronics circuit simulator software released under GPL license. 
It gives the possibility to set up a circuit with a graphical user interface built upon QT libraries 4 and 



1 


http://qucs.sf.net 






2 


http://www.gnu.org/licenses/gpl.html 




3 


https://github.com/zonca/python-qucs 




4 


http://trolltech.com/products 
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simulate the signal and noise behaviour of the circuit. It is an alternative to well known Berkeley 
SPICE 5 and Agilent ADS. 

QUCS is coded in C++ and uses extensive class inference in order to facilitate the implemen- 
tation of new components. The main QUCS modules are: 

qucsator the command line circuit simulator, which reads a circuit description in a predefined 
ASCII format, called netlist, and outputs an ASCII format results file. 

GUI Graphical User Interface, which is completely independent and makes it possible to draw a 
circuit using a library of devices or file defined components. The GUI automatically builds 
the netlist, runs qucsator, parses the results file, and allows the user to easily produce 
tables and plots. 

2.1.1 Solving method for AC simulations 

QUCS simulator can be used for simulating electrical circuits in time or frequency domain, for 
linear or non linear analysis. The simulation used in this article is based on an AC analysis of the 
circuit. The AC analysis is done by using a Modified Nodal Analysis (MNA) method, see QTIf ]. 
MNA is based on solving the circuit equations in the form of a linear system: 



The A matrix is only dependent on characteristics of the components, the x vector contains 
the unknown quantities (node voltage and current through independent voltage sources), and z 
represents the vector of sources (voltage and current sources). Simulation is done simply by solving 
the system. 

The main strength of the MNA method is the stamp approach: each component is characterized 
by four small matrices and one vector that are directly inserted into the A and z matrices. Thus the 
matrices are built from small building blocks without general knowledge of the full circuit topology. 

2.2 Bandpass characterization 

Radiometric devices are receivers capable of measuring the power of tiny microwave signals by 
providing strong amplification, typical gains of about 10 7 in terms of power, and detection through 
square law detector diodes. A wide variety of receivers have been used in microwave and millime- 
ter wave astrophysics, from simple total power systems, to Dicke-switched receivers, to various 
correlation schemes . The simplest radiometer, Fig. [I], is composed by a feed horn which con- 
veys the signal to an amplification stage, a bandpass filter which defines the frequency range and a 
diode which outputs a DC voltage that is proportional to the power of the incident electromagnetic 
signal and can be easily measured and digitized. 

Their bandpass is defined as the normalized gain as a function of frequency, see 0] : 



[A]x = z. 



(2.1) 




(2.2) 



where: 



5 http://bwrc.eecs.berkeley.edu/Classes/IcBook/SPICE/ 
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• V out is the radiometer voltage output 

• G is the calibration factor in units [V/W] 

• g(v) is the normalized bandpass 

• Pi„ is the radiometric input power 

• V noise is the voltage due to the radiometer noise 

In radiometer frequency response modeling we can devise two main cases: characterization of 
only the RF section of a radiometer or end-to-end including the diode. 

In the first case QUCS just combines the S-parameters of each device taking into account 
multiple reflection and outputs the S-parameters of the assembly, S21 is the bandpass. 

In the second case the diode is quadratic so that the transfer function is linear between in- 
put power and output volts, so that it is necessary to run AC simulations at two different input 
temperatures and compute the gradient as a function of frequency: 

Vout 2 {v) -Vouhiv) 
{) Win 2 (v)-Wm(v)' K ' 

Each component can be modeled by measured data or simulated; measurements data consists 
of a Touchstone 6 file with the four S-parameters and an optional noise figure, both as a function of 
frequency. 

In case it is necessary to estimate the bandpass response on a larger frequency span than 
effectively measured, it is easier to extrapolate on the model than an end-to-end measurement, 
since the low and high frequency cutoff sections of the response on the measurement are a non- 
trivial combination of each device frequency response. An example is given in section ^3] for the 
extrapolation of Planck-LFI OrthoMode Transducers response. 

2.3 Rectangular waveguide model 

Since QUCS is Free Software, it is possible to access the source code and implement functional- 
ities which are missing. In order to model LFI rectangular waveguides we have implemented a 
rectangular waveguide model which was then submitted to the QUCS repository and merged into 
the following QUCS release. 

In this section we describe the implementation of the waveguide device in QUCS. 

2.3.1 Implementation in QUCS 

The rectangular waveguide model implementation is based on the S-parameter components and has 
been quite straightforward thanks to the modular conception of QUCS. The code consists of about 
600 lines including extensive comments and the GUI. 

The waveguide component was validated using Agilent ADS results with a reference waveg- 
uide model defined as follows: 

• single rectangular waveguide component 

6 Touchstone is a standard text format for describing S-parameters and noise properties as a function of frequency 



-5- 




• WR-22 (same as 44 GHz LFI waveguides) 

• 800 mm length 

• material is gold plated stainless steel 

• simulation temperature is 20 °C 

Both S parameter and noise simulation voltage results were compared and the matching is at 
the level of a tenth of a dB, figure @ shows S21 as a function of input frequency of ADS and QUCS 
over-plotted. 

We also added the embedded computation of the resistivity of gold, stainless steel and alu- 
minum based on device temperature. This feature was implemented using different empirical for- 
mulae in different temperature ranges. 



2.4 Batch processing 

In the context of LFI bandpass modeling we developed a set of Python tools for automating the 
procedure of preparing, running and exporting data from the simulations. 

They are based on the software package SciPy 7 , a general purpose scientific and numerical 
software package for Python, in particular on: 

• NumPy: multidimensional array package 

• Matplotlib: 2D plotting library 

• Ipython: interactive Python console 

7 http://www.scipy.org 
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The programs are implemented following Object Oriented Programming and are strictly inde- 
pendent, each stage runs separately by reading inputs and writing results in standard ASCII files, 
the stages are: 

XLS2S2P: Format conversion from data stored in an Excel (or Gnumeric) spreadsheet, easier to 
build and maintain, to Touchstone, used by QUCS. 

python-qucs: Batch simulation run using qucsator (QUCS GUI is never launched) and 
exporting of results from QUCS to ASCII format; a simulation of the 44 GHz LFI channels 
takes about 5 minutes on a common laptop. 

ba_lib: Bandpass analysis tool, for batch processing and interactive analysis: 

• gain bandpass from simulations at 2 different temperatures 

• normalization 

• plotting 

• comparison with swept source measurements and ADS simulation 

• export to text format 

These stages were built for PLANCK-LFI simulations but implement several general utilities 
that can be assembled to build a simulation pipeline for other experiments. 

2.4.1 python-qucs 

python-qucs, the most general tool we developed, was released under GPL on Github 8 . It was 
designed to easily iterate between testing a circuit on the QUCS GUI and running batch simulations, 
so that a system can be debugged interactively at several development stages. The simulate 
module parses the netlist produced interactively by the QUCS GUI and runs user defined regular 
expressions to create and save to disk the target set of netlists, then runs qucsator on each 
of them. The extract and plot modules allow then to read the custom QUCS file format, 
manipulate and plot the results. 

3. Application examples: COFE and Planck-LFI 

In this section we present the application of the QUCS bandpass model to bandpass modeling of 
two radiometric instruments for CMB measurements: COFE (COsmic Foreground Explorer) and 
PLANCK-LFI (Low Frequency Instrument). 

3.1 COsmic Foreground Explorer 

COFE, [ ]I3T| , is a balloon-borne microwave polarimeter under construction at the University of Cal- 
ifornia, Santa Barbara and designed to measure low-frequency dominant diffuse polarized galactic 

8 https ://github. com/zone a/py thon-quc s 
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foregrounds. The first 24-hour flight, foreseen in summer 2011 from Fort Sumner, New Mex- 
ico, will produce temperature and polarization maps at 10 and 15 GHz of 59% sky with median 
aggregate sensitivity of 92 [iK/deg 1 for temperature and 77 /j.K/deg 2 for polarization. 

COFE is comprised of an off axis Gregorian telescope with a reflecting quarter wave plate 
polarization modulator at the co-focal point of the two focusing optical elements. Using an external 
chop to overcome 1/f noise allows COFE to use a total power radiometer. Each receiver chain has 
the same design and consists of a cryogenically cooled (~ 20 K) InP MMIC LNA 9 and provides 
roughly 30 dB of gain with a noise temperature of 8K. Following the LNA in the cryostat are 2 
commercially available GaAs 10 MMICs, a bandpass filter and a square law detector diode. The 
LNA is connected to the ambient temperature components via low loss coax ( IdB). 



Cryo Amp 


Warm Amplifiers Filter 






-40 dB Jj ^ -3dB 


-IdB _?~ m ± r -6dB ^ 





Figure 3. COFE 10 GHz radiometer VNA test setup for bandpass characterization 



In order to estimate the bandpass response of the RF part of the 10 GHz receiver we performed 
a dedicated warm Vectorial Network Analyzer (VNA), Rohde&Schwarz ZVA40, test on a setup 
composed of, see figure 

• 40 dB attenuation at the input, necessary to reduce the input power generated by the VNA to 
the expected power level in flight 

• cryogenic low noise amplifier, ~ 30 dB gain 

• 3 dB attenuator, which is the expected loss of the coaxial cable that will connect the cryogenic 
section to the warm section 

• first warm amplifier, ~ 20 dB gain 

• 1 dB attenuator, which will be present also in flight, useful to provide isolation to the ampli- 
fiers 

• second warm amplifier, ~ 20 dB gain 

• 6 dB attenuator designed to provide isolation between the bandpass filter and the amplifying 
stage 

• bandpass filter, designed to provide 4 GHz bandwidth and strong out of band rejection 

9 Indium Phosphide Monolithic Microwave Integrated Circuit Low Noise Amplifier 
"'Gallium Arsenide 
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All of the components were tested individually and as an assembly with and without the filter 
on the VNA measuring all S-parameteres. 

We built a COFE 10 GHz radiometer model in QUCS using the measured S-parameters of each 
component and ran simulations in order to have an independent estimate of the bandpass response. 




Figure 4. COFE 10 GHz radiometer test and model results without bandpass filter, this includes the 4(MB 
attenuation as part of the setup, the true amplification of the radiometer is almost 70 dB 

Figures |] and |5J show the comparison of modeled and measured bandpasses, without and with 
the filter, respectively. The agreement within the band is better than 1.5 dB for the radiometer with 
the filter, which is at the same level of the expected accuracy of the measurement. 

The bandpass model proved to give a reliable estimate of the bandpass response of the COFE 
10 GHz radiometer. It was also applied to support the filter design and the estimate of the out-of- 
band rejection. 

In the future we plan to use the COFE bandpass model and python-qucs for automating 
the process of estimating the best combination of the available components. The COFE 10 GHz 
radiometer has only 3 receivers, but for each of them the cryogenic amplifier can be coupled to any 
2 of the 6 available warm amplifiers, their order matters, so the total number of combinations is 90. 
We will estimate the bandpass using the model for all the combinations and select just the best in 
terms of bandwidth for performing hardware testing. 

3.2 Planck-LFI 

As a second example, we consider the application of our model to the simulation of the bandpass 
response of the PLANCK-LFI receivers. We provide an overview of the LFI radiometer scheme, 
present an example of extrapolation of the response of a device from limited measured data and 
then compare results of simulations run with QUCS and Agilent ADS [0. 
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COFE Bandpass response with filter 

40 1 1 1 1 1 




Frequency [GHz] 

Figure 5. COFE 10 GHz radiometer test and model results with bandpass filter 

LFI jni , ^] is a radiometer array mounted in the focal plane of PLANCK satellite telescope, 
receiving microwave photons at 30, 44 and 70 GHz. 

The complete LFI is an array of 1 1 Radiometer Chain assemblies (RCA), 6 with a central 
frequency of 70 GHz, 3 of 44 GHz and 2 of 30 GHz. Each RCA is composed of (Fig. §): 

• a Feed Horn looking at the Sky signal 

• an OrthoMode Transducer (polarisation splitter) 

• two reference horns optically coupled to two reference loads at about 4K 

• a cold (~20K) pseudo correlation stage (Front End Module) providing ~ 35 dB amplification 

• four waveguides connecting the cold module to the warm backend 

• a second RF amplification stage and bandpass filter 

• a square law detector diode 

• post-detection DC electronics stage at 300 K (Back End Module) 
3.3 Component data extrapolation 

LFI OrthoMode Transducers (OMT), see QT3|], are passive components based on waveguides that 
split the incoming radiation into its orthogonal polarized components. Because the OMT design 
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Figure 6. Picture of a 30 GHz Radiometer Chain Assembly (RCA) mounted for testing in Milano 



was exactly the same for 30 and 44 GHz units, it is possible to scale at 30 GHz the measurements 
performed at 44 GHz, and vice versa. 

In particular we estimated the band between 21.3 GHz and 26.5 Ghz (which was outside the 
limit for direct measurements) by rescaling 44 GHz measurements which were performed on a 
larger band, from 33 to 50 GHz. 

The extrapolation (see figure [7|) was done after normalizing the frequency axis on the central 
frequency, so that all the curves have a central frequency of 1. 

3.4 Model and results 

Implementation, excluding the waveguides, faithfully follows the ADS model already described in 
[f|] and [01. For reference we show the model schematic in figure || 

In figure we show the comparison of estimated bandpasses by ADS and QUCS. As expected 
QUCS and ADS results are identical for measured components and differ only for the waveguide 
simulation, the difference is about ±2 dB, which is at the same level of the expected bandpasses 
accuracy. 

4. Conclusion 

In this paper we have presented a method and the software tools capable of estimating radiometer 
bandpasses from individual component measurements and analytical simulations based on QUCS 
and Python. The approach is general and can be used during the design phase of an instrument, 
gradually replacing analytical simulations, based on requirements, with real data as soon as they 
become available, in order to continously support frequency response estimation with a software 
simulation for comparison and validation. 

We presented two application scenarios: in the COFE experiment the model was successfully 
validated by comparison with VNA measurements and is going to be used to optimize component 
matching; in the PLANCK-LFI mission the model replaced a model based on ADS thanks to the 
possibility of easily batch processing simulations with Python. 




Figure 7. Extrapolation of the S21 30 GHz measurements compared with 44 GHz measurements normalized 
to their own central frequency. The extrapolation was necessary between 21.3 GHz and 26.5 GHz (black 
vertical line). At high frequency the cutoff was not characterized because the BEM bandpass filter has strong 
rejection in that region 
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Figure 8. Model schematic: each of the 44 LFI channels is modeled independently as the above model. 
Each component is represented by its measured behavior as a function of frequency. Waveguides instead are 
simulated analytically given their dimensions, length and temperature gradient along their length. 

python-qucs is a flexible and easy-to-use software tool for making QUCS modeling scale 
with the number of channels required by the next generation CMB experiments. It provides both 
powerful batch simulation capabilities for running many channels in different conditions and an 
interactive debug session for better insight on any issue. 

Therefore, QUCS frequency response modeling and python-qucs offer both performance 
estimation and measurements support in the context of radiometers bandpass estimation, with low 
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RCA27 channel S-ll 




Frequency[GHz] 

Figure 9. Comparison between ADS and QUCS bandpasses on LFI27S-1 1 at 30 GHz 

cost in terms of budget, since both pieces of software are open-source, and in terms of time, since, 
usually, individual component measurements are already available and the tools strongly reduce 
the time needed to build and maintain a model. 
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